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ABSTRACT 

Aims. The focusing performance of X-ray optics (conveniently expressed in terms of HEW, Half Energy Width) strongly depend on 
both mirrors deformations and photon scattering caused by the microroughness of reflecting surfaces. In particular, the contribution 
of X-ray Scattering (XRS) to the HEW of the optic is usually an increasing function H(E) of the photon energy E. Therefore, in 
future hard X-ray imaging telescopes of the future (SIMBOL-X, NeXT, Constellation-X, XEUS), the X-ray scattering could be the 
dominant problem since they will operate also in the hard X-ray band (i.e. beyond 10 keV). In order to ensure the imaging quality at 
all energies, clear requirements have to be established in terms of reflecting surfaces microroughness. 

Methods. Several methods were proposed in the past years to estimate the scattering contribution to the HEW, dealing with the surface 
microroughness expressed in terms of its Power Spectral Density (PSD), on the basis of the well-established theory of X-ray scattering 
from rough surfaces. We faced that problem on the basis on the same theory, but we tried a new approach: the direct, analytical 
translation of a given surface roughness PSD into a H(E) trend, and - vice versa - the direct translation of a H(E) requirement into a 
surface PSD. This PSD represents the maximum tolerable microroughness level in order to meet the H(E) requirement in the energy 
band of a given X-ray telescope. 

Results. We have thereby found a new, analytical and widely applicable formalism to compute the XRS contribution to the HEW 
from the surface PSD, provided that the PSD had been measured in a wide range of spatial frequencies. The inverse problem was also 
solved, allowing the immediate evaluation of the mirror surface PSD from a measured function H(E). The same formalism allows 
establishing the maximum allowed PSD of the mirror in order to fulfill a given H{E) requirement. Practical equations are firstly 
developed for the case of a single-reflection optic with a single-layer reflective coating, and then extended to an optical system with 
A' identical reflections. The results are approximately valid also for multilayer-coated mirrors to be adopted in hard X-rays. These 
results will be extremely useful in order to establish the surface finishing requirements for the optics of future X-ray telescopes. 
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1. Introduction 

The adoption of grazing-incidence optics in X-ray telescopes in 
the late 70s allowed a great leap forward in X-ray astronomy 
because they endowed the X-ray instrumentation with imag- 
ing capabilities in the soft X-ray band {E < 10 keV). The 
excellent performances of the soft X-ray telescopes ROSAT 
(Aschenbach [T988] l. Chand ra (Weisskopf |2003] ) and Newton- 
XMM (Gondoin et al. 1998) are well known. 

To date, the utilized technique to focus soft X-rays consists 
in systems of double-reflection mirrors with a single layer coat- 
ing (Au, Ir) in total external reflection at shallow grazing inci- 
dence angles. In this case, the incidence angle 0i (as measured 
from the miiTor surface) cannot exceed the critical angle for total 
reflection, otherwise the mirror reflectivity would be very low. 
The critical angle is inversely proportional to E, the energy of 
the photons to be focused. Using Au coatings, for instance, the 
incidence angle cannot exceed ~ 0.4 deg for photon energies 
E ^ lOkeV. 

An extension of this technique to the hard X-ray energy band 
(£■ > 10 keV) can be pursued by combining long focal lengths 
(> 10 m), very small incidence angles (0.1 H- 0.25 deg), and 
wideband multilayer coatings to enhance the reflectance of the 
mirrors at high energies (Joensen et al. |1995t Tawara et al. |1998l ). 
A very long focal length is hardly managed using a single 
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spacecraft, therefore the optics and the focal plane instruments 
should be carried by two separate spacecrafts in formation- 
flight configuration. This is the baseline for the future X-ray 
telescopes SIMBOL-X (Pareschi & Ferrando l2006b and XEUS 
(Parmar et al. |200 4). Other hard X-ray imaging telescopes of 
the future are NeXT (Ogasaka et al. 2006) and Constellation-X 
(Petre et al. l2006] l. 

The focusing and reflection efficiency of X-ray optics can be 
tested and calibrated on ground by means of full-illumination 
X-ray facilitie s hke PANTER (Briiuninger et al. l2004l 
Freyberg et al. 120061 ). successfully utilized in the last years to 
calibrate the optics of a number of soft X-ray telescopes. The 
PANTER X-ray facility now allows testing in soft (0.2 H- 10 keV) 
and hard (15 H- 50 keV) X-rays multilayer-coated optics pro- 
totypes for future X-ray telescopes (Pareschi et al. 120051 
Romaine et al. I2005I I. The source distance finiteness causes 
some departures of the optic performances, with respect to 
the case with the source at astronomical distance: elTective 
area loss, different incidence angles on paraboloid and hyper- 
boloid, focal length displacement, a slight focal spot bluiTing 
(Van Speybroeck & Chase 1972). However, there effects can 
be quantified and subtracted from experimental data. After this 
treatment, the focusing-concentration performances of the optic 
can be experimentally characterized as a function of the inci- 
dent photon energy, in terms of Half-Energy Width (HEW) and 
Effective Area (EA). 
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The focusing performance, in particular, is altered by mirror 
deformations that may arise in the manufacturing, handling, in- 
tegration, positioning processes. The consequent imaging degra- 
dation can be calculated from the measured departures of the 
mirrors from the nominal profile, by means of a ray-tracing pro- 
gram. As long as the geometrical optics approximation can be 
applied, the effect is independent of the photon energy. The fig- 
ure errors contribution to the HEW can be also directly measured 
using a highly coUimated beam of visible/UV light in a precision 
optical bench. In this case, however, the light diffraction has to 
be carefully estimated and subtracted. 

Another drawback is the X-ray scattering (XRS) caused by 
the microroughness of reflecting surfaces (Church et al. 119791 
Stearns et al. 1998; Stover 1995; and many others). The XRS 
spreads a variable fraction of the reflected beam intensity in the 
surrounding directions: the result is the effective area loss in the 
specular direction (i.e. in the focus) and a degradation of the 
imaging quality. The XRS is an increasing function of the pho- 
ton energy; due to the impact that the XRS can have on astro- 
nomical X-ray images quality, the height fluctuations rms of the 
mirror surface should not exceed few angstroms. Loss of effec- 
tive area is also caused by interdiffusion of layers in multilayer 
coatings, which enhances the X-ray transmission and absorption 
throughout the stack. On the other hand, an uniform interdiffu- 
sion does not cause X-ray scattering (Spiller [T994| l. hence it does 
not contribute to the focusing degradation. 

The microroughness of an X-ray mirror can be measured on 
selected samples using several metrological instruments, each of 
them sensitive to a definite interval of spatial scales i. Long Trace 
Profilometers (10 cm > / > 0.5 mm: Takacs et al. 1999), opti- 
cal interference profilometers (5 mm > / > 10 fim) and Atomic 
Force Microscopes (100 /im > I > 5 nm) can be suitable in- 
struments to provide a detailed profile characterization of X-ray 
mirrors surface. It is convenient to present the deviation of sur- 
face from the ideality in terms of Power Spectral Density (PSD), 
because its values do not depend on the measurement technique 
in use (see llSO lOIlO Standardl l. In addition, the XRS diagram, 
and consequently the HEW, can be immediately computed from 
the PSD at any photon energy (Church et al. |1979l l. 

In the past years, several approaches were elaborated to 
relate a mirror PSF (Point Spread Function) to the PSD 
of its surface. Among a wealth of works, we can cite 
(De Korte et al. |1981| l the assumption of a Lorentzian model 
for the PSD to fit the mirror PSFs at some photon energies, 
allowing the derivation of two parameters (roughness rms and 
correlation length) of the model PSD. Christensen et al. ( I1988I I 
perform a fit of experimental high-resolution XRS data dealing 
with the surface correlation function. Harvey et al. (1988) re- 
late the PSF of Wolter-I optics to the parameters of an exponen- 
tial self-coiTelation function along with a transfer function-based 
approach. Willingale (,1988) derived the surface PSD of a miiTor 
from the wings of a few PSFs, measured at PANTER at some 
soft X-ray photon energies. O'Dell et al. ( 1993) interpret the PSF 
of a focusing mirror on the basis of surface roughness and partic- 
ulate contamination. Zhao & Van Speybroeck (I2003I I construct 
from the PSD of a focusing mirror a model surface and compute 
the X-ray scattering PSF from the Fraunhofer diffraction theory. 

In the present work that problem is faced in a new and differ- 
ent way, looking for a general and simple link between measured 
roughness and mirror HEW. More precisely, we considered the 
following question: for an X-ray grazing-incidence optic, what 
is the maximum acceptable PSD of the surface that fulfills the 



angular resolution (HEW) requirements of the telescope, in all 
the energy band of sensitivity! 

In this work we shall give a definite answer to this question. 
In the sect. |2] we shall summarize the causes of imaging degra- 
dation. In the sect. [3] we show how to evaluate H(E), the XRS 
contribution to the HEW of a focusing mirror at the photon en- 
ergy E, from any surface microroughness PSD, measured over a 
very wide range of spatial frequencies. We shall see in the sect.|4] 
that for the special class of fractal surfaces we can even relate 
the power-law indexes of PSD and HEW, and in the sect. |5] we 
see how to treat the other cases. Then we prove in the sect.|6]that 
the formalism can be reversed, providing thereby an independent 
evaluation of the surface PSD from an analytical calculation over 
H(E), and in the sect. |7] we extend the results to focusing mirrors 
with more than one reflection. Finally, an example of computa- 
tion is provided in the sect. [8] 

2. Contributions to the imaging degradation 

We shall henceforth indicate with A the wavelength of photons 
impinging on the mirror, and we shall consider the HEW as a 
function of A instead of the photon energy E. For isotropical re- 
flecting surfaces in grazing incidence, the X-ray scattering dis- 
tribution lies essentially in the incidence plane, so we denote the 
incidence angle on the mirror as 0i and the scattering angle as 6^, 
both measured from the surface plane (a schematic of the scatter- 
ing geometry is drawn in fig. [TJ. If we do not consider the optic 
roundness errors, the longitudinal deviations from the nominal 
profile of a focusing mirror can be classified on the basis of their 
typical length /. According to De Korte et al. ( |1981| l. they are: 

1 . Power errors: errors with / equal to the mirror length L. They 
consist in a single-concavity deformation of the profile with 
respect to the nominal one. 

2. Regularity errors: errors in the spatial range from 0.1 L < 
/ < 0.5 L. 

3. Surface roughness: surface defects with I <Q.l L. 

However, other criteria were also formulated to separate fig- 
ure errors from roughness. Consider a single Fourier component 
of the surface profile with wavelength / and root mean square cr. 
That Fourier component is dominated by figure error if it fulfills 
the condition (Aschenbach l2005t 

47rsin6iiO- > /i. (1) 

Otherwise, it is dominated by microroughness. In other words, 
surface defects within the smooth-surface approximation can be 
mainly considered as microroughness. To understand the impor- 
tance of this approximation, we write the optical path difference 
As of X-rays reflected by two points of the surface with a hori- 
zontal spacing I and vertical spacing & - 2 V2cr (for optically- 
polished surfaces, cr is a increasing function of I, and usually 
o- « 10"^^/) as 

As - [(cos 0s - cos 6i) + (T(sin 0i + sin 0s) (2) 

that, for small incidence angles, becomes 

As = fsin 0i(0, - 0i) + d-(0i + 0,). (3) 

If that component is responsible for X-ray scattering, it has 
to be As X A, to cause the diffraction from surface features with 
a / spacing and & height. Conversely, the "figure errors", which 
are treated with the methods of the geometrical optics, should 
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r^iectedX-ra^s 

0. = 6^ 



Fig. 1. The geometry of X-ray scattering: the strictly speaking 
"reflected" rays (i.e. in the focus direction) are characterized by 
the equality 0, = 0,, the others are scattered apart. The rough sur- 
face is a simulated one, assuming a PSD with power-law index 
n - 2.4 (see sect.|4|i. 



be characterized by the inequality hs » A. Note that this condi- 
tion becomes similar to the eq.[T]in the limit \6\ - — > 0. The 
application of this criterion and of the subsequent X-ray scatter- 
ing theory requires the incident radiation to be spatially coher- 
ent over the spatial scale /, so that the properties of the reflected 
wavefront are determined only by the coherence properties of 
the mirror surface. This in turn requires the angular diameter of 
the source 0s to fulfill the inequality (Holy et al. II 9991) 



/ sin ft 



(4) 



This equation sets a maximum to the values of I that can be used 
in the application of the results presented in this work. The limi- 
tation can affect X-ray sources at finite distance, like those used 
for X-ray optics calibrations in full-illumination setup. For very 
distant astronomical X-ray sources, the conditionals met even 
for larger /, up to Z a; L. 

It is worth pointing out that, for a given reflecting surface, 
the separation of figure errors from microroughness is strongly 
affected by the incidence/scattering angles. In fact, even for large 
I, Ai can become comparable with /I, if Q\ and are sufficiently 
small: thus, the spatial wavelength window of interest for X-ray 
scattering can shift to the large / domain (or, equivalently, to the 
range of low spatial frequencies / = 1//), provided that the con- 
dition|4]is fulfiUed. 

Let us now consider how to separate the figure and scattering 
terms in HEW data. In absence of XRS, the mirror PSF would 
be independent of the energy and due only to figure errors (i.e. 
in the approximation of the geometrical optics). The resulting 
HEW would be also constant. Instead, due to the XRS, the figure 
PSF is convolved with the X-ray scattering PSF to return the 
PSF(^) being measured (Willingale [T9881 Stearns et al. [19981 
and many others). 



(5) 



The resulting HEW will depend on the photon wavelength, 
as it does the PSF. In order to isolate the scattering term from the 
total PSF a deconvolution should be carried out, provided that 
the PSFfig is known. However, if we assume that the XRS and 
the mirror deformations are statistically independent, the total 



HEW can be approximately calculated as the squared sum of the 
two contributions: 



HEW^(A) 



HEWl„ + H^iX). 



(6) 



An estimation of HEWf,;, can be obtained: 



1. from the application of a ray-tracing code to several mea- 
surements of the mirror profile, 

2. from reliable extrapolation of the HEW{A) curve to £ — > 0, 
in absence of low-energy diffraction effects like dust contam- 
ination, studied in detail by O'Defl et al. (119931) 

3. from a direct measurement of the HEW in visible/UV light, 
provided that the diffraction at the mirror edges can be reli- 
ably calculated and subtracted. 

Once known the measured HEW(/l) experimental trend and 
the HEWfig term, the eq. |6]can be used to isolate the scattering 
contribution from the experimental HEW trend: we shall prove 
in the next section that the H{A) function is immediately related 
to the reflecting surface ID Power Spectral Density (PSD) P{f) 



p(f> = i 



Jo 



z(x)e 



-2mf 



dx 



(7) 



where z(x) is a height profile (of length L) of the miiTor, mea- 
sured in any direction (Stover [T9951 l: the surface is assumed to 
be isotropic, and the spectral properties of the profile to be rep- 
resentative of the whole surface. The PSD is often measured in 
nm^ units, and for optically-polished surfaces it is usually a de- 
creasing function of the frequency /. 

PSD measurements have always a finite extent [/min,/max], 
determined by the length and the spatial resolution of the mea- 
sured profile. As well known, the surface rms <t is simply com- 
puted from the PSD by integration over the spatial frequencies /: 



P(f)df 



note that the integration range should always be specified. 



(8) 



3. Estimation of H{A) for single-reflection focusing 
mirrors 

3.1. Single-layer coatings 

Firstly, we suppose the mirror to be plane and single-layer 
coated. For a surface with roughness rms o", the specular beam 
intensity obeys the wefl-known Debye-Waller formula 



R - Rp exp 



leTrV^sin^^i 
^2 



(9) 



here Rf is the reflectivity at the grazing incidence angle fti, as 
calculated from Fresnel's equations (zero roughness). However, 
it should be noted in the eq. |9]that neither the spatial frequen- 
cies range where the PSD should be integrated is specified, nor 
the separation between reflected and scattered ray is clearly in- 
dicated: these ambiguities can be solved as follows. 

Let us derive the total scattered intensity Is from the con- 
servation of the energy: for smooth surfaces, i.e. fulfilling the 
inequality 2cr sin 6^ <K A, we can approximate 



loRf 



1 - exp 



167rV2sin^6»i 
l2 



leTrVsin^^i 
IoRf j2 ^-(10) 
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In grazing incidence. X-ray scattering lies mainly in the inci- 
dence plane. Moreover, the normalized scattered intensity per ra- 
dian at the scattering angle (either 6*5 > 6i or < Oi) is related 
to the PSD along with the well-known formula at first-order ap- 
proximation (Church et al. 11979', Church & Takacs 1986 ), valid 
for smooth, isotropic surfaces and for scattering directions close 
to the specular ray (i. e. \0s - Oi\ <K 9i), 



(11) 



1 d/s 16n^ , 

where P(f) is the Power Spectral Density of the surface (eq. |7]l 
and Iq is the flux intensity of the incident X-rays. If the scattered 
intensity is evaluated at the scattering angle 9s, the PSD can be 
immediately evaluated as a function of the spatial frequency /: 



J cos 9i - cos 0s 



sin0i(0s-0i) 



(12) 



A A 

In the eq. [T2]the approximation was justified by the assumption 
l^s - 9i\ <K 9\ and the negative frequencies are conventionally 
assumed to scatter at 9^ < 9i. the assumed approximations make 
the XRS diagram symmetric, because the PSD is an even func- 
tion. 

For a single-reflection mirror shell, the extension of the for- 
mulae above-mentioned is straightforward by regarding \9s - 9\\ 
as the angular distance at which the PSF is evaluated. The fo- 
cal image is the superposition of many identical XRS diagrams 
on the image plane, generated by every meridional section of the 
mirror shell: since a n angle rotation of every meridional plane of 
the shell sweeps the whole image plane, the scattered intensity is 
spread over a n angle. The integration on circular coronae used 
to compute the mirror PSF (at positive angles) compensates this 
factor multiplying the XRS diagram by 2n (De Korte et al. 198T]|. 
The remaining 2-fold factor accounts for the negative frequen- 
cies in the surface PSD. We shall henceforth suppose that the fac- 
tor 2 is embedded in the PSD definition. Therefore, the eqs.fTTI 
and[12]can be used to describe the XRS contribution to the PSF. 

We are now interested in the scattered power at angles larger 
than a definite angle a measured from the focus. Due to the steep 
fall of scattering intensity for increasing angles, the integral has 
a finite value 



I[\9s-9i\>a]^ -^d9s. 

JSi+a It's 

Combining eqs.fTTI and [T3l one obtains: 



I[\9s-9i\>a]^IoRF- 



16n^ sin^ 9i 
l3 



f P(f)d9s 



(13) 



(14) 



with respect to the definition used in the eqs. |7]and[TT] a factor 
2 was included in the PSD. The upper integration limit corre- 
sponds to a photon back-scattering: at first glance, this seems 
to violate our small-scattering angle assumption (eqs. [TT| and 
[T2I 1. but it should be remembered that only the angles close to 9i 
contribute significantly to the integral in eq. [13] hence its value 
should not be significantly affected by a particular choice of the 
upper integration limit. After a variable change from 9s to / 
(eq. [T2I1. the eq. [T4l becomes (approximating cos^i ^ 1 in the 
upper integration limit): 



I[\9s-9i\>a]^IoRF 



16n^ sin^^i 
^2 



P(f)df 



(15) 



where /o - asm9i/A is the spatial frequency corresponding to 
the scattering at the angle a. As expected, this equation equals 



the integrated scattering according to the eq. [TOl provided that 
we identify / [\0s - 9i\ > a] with Is, and the squared roughness 
rms with 



P(f)df. 



(16) 



The eq.[T6]is in agreement with the eq.[8] but it states clearly the 
window of spatial frequencies involved in the XRS. Therefore, 
for a definite angular limit a the "reflected beam" intensity can 
be simply calculated by using the Debye- Waller formula, pro- 
vided that cr^ is computed from the PSD integration beyond the 
frequency /o, which corresponds to an X-ray scattering at a. The 
upper integration limit is a very high frequency (close to 1/A): 
hence, the atomic structure of the surface is not important in the 
integral of the eq. [16] Moreover, considering that the PSD trend 
for optically-polished surfaces decreases steeply for increasing 
/, the largest contribution to the integral should be given by the 
frequencies close to /q. 

Now we can evaluate H(A), the scattering term of the HEW. 
For simplicity, in the following we will suppose that the HEW 
is obtained from the collection of all the reflected/scattered pho- 
tons: this allows us to avoid problems related to the finite size of 
the detector, and to extend the surface roughness PSD up to very 
large spatial frequencies. By definition, H{A) is twice the angu- 
lar distance from focus at which the integrated scattered power 
halves the total reflected intensity: 



I[\9s-9,\>a] = -IqRp 

we immediately derive, from the eq.|9] 



exp - 



16n^o-^ sin"- 9i 



(17) 



(18) 



where cr^ has now the meaning as per the eq. [16] Solving the 
eq. [T8]for cr^ and equating to the integral of the PSD, 



P(f)df 



A^ In 2 
167r2 sin^ 9i ' 



(19) 



once known the PSD from topography measurements over a 
wide range of spatial frequencies, the PSD numerical integra- 
tion in the eq. [19] allows to recover /q. In turn, fit is related to 
H(A) through the eq.[T2] that we write in the following form 



H(A) 



2Afo 
sin 9i ' 



(20) 



where H is measured in radians. Note that the condition 
H(X) <K 9i is very important, for the eq. [20] to hold. Small 
scattering angles and grazing incidence are also very important 
for the considerations that follow. 



3.2. Multilayer coatings 

The obtained result (eq. [T9b can be extended to mirrors with 
multilayer coatings, used to enhance the grazing incidence re- 
flectivity of mirrors in hard X-rays {E > 10 keV). In general, 
the multilayer cannot be characterized by means of a single 
PSD, due to the evolution of the roughness throughout the stack 
(Spiller et al. 1993; Stearns et al. 1998). Moreover, due to the 
interference of scattered waves at each multilayer interface, the 
final scattering pattern is more structured than eq.[TT] with peaks 
whose height depends on the phase coherence of the interfaces 
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— p(n), dependence on Kp 
y(n), dependence on X 




2.0 2.5 
PSD power-law index 



3.5 



Fig. 2. Dependence of the spectral exponents for different in- 
dexes « of a power-law PSD, for a single-reflection focusing 
mirror. In the forbidden region (n > 3) y would be negative. 



(Kozhevnikov '2003 ") . The HEW term can be computed numeri- 
cally from the XRS diagram. 

In order to extend the eq. [T9]to mirrors coated with a graded 
multilayer, we have to assume the additional requirements: 

1. the PSD is constant and completely coherent throughout 
the multilayer stack: i.e., the deposition process does not 
cause additional roughness and replicates simply the profile 
of the substrate. Therefore, all the PSDs and all the cross- 
correlation between interface profiles equal the PSD mea- 
sured at the multilayer surface. This is often observed in 
the I > lOjum regime (Canestrari et al. I2006I I. where most 
of frequencies /o fall when the incidence angle is less than 
0.5 deg. Most of microroughness growth, indeed, takes place 
for 10//m> I > 0.1 yum. 

2. the multilayer reflectivity RAiO;) at the photon wavelength A 
changes gradually over angular scales of H{A). Ideally, this 
condition should be fulfilled by wideband multilayer coat- 
ings for astronomical X-ray mirrors. 

Under these hypotheses, a quite tedious calculation reported 
in appendix lAl shows that the eq. [T9]can be approximately ap- 
plied also with multilayer coatings. The following developments 
also apply in that case. 



4. H{A) for a fractal surface 

We apply now the equations [19] and |20] to the typical (monodi- 
mensional) PSD model for optically-polished surfaces, a power- 
law (Church il988i l 



K„ 



(21) 



where the power-law index « is a real number in the interval 
1 < « < 3 and K„ is a normalization factor. A power-law PSD is 
typical of a. fractal surface, and it represents the high-frequency 
regime of a K-correlation model PSD (Stover [T995l l. This model 
exhibits a saturation for / — > that avoids the PSD divergence. 
In practice, the fractal behavior dominates in almost all spatial 
frequencies of interest for X-ray optics. 

There are interesting reasons for which n can take val- 
ues on the interval (1:3). In fact, for a surface in the 3D 
space, n is related to its Hausdorff-Besicovitch dimension D 
(Barabasi & Stanley [1995]l along with the equation n - 1 - 2D 



(see Church [19881 Gouvet [T996] l. The restriction 1 < « < 3 for 
a fractal surface is therefore necessary to have 3 > D > 2. 

A power-law PSD is particularly interesting because the in- 
tegral on left-hand side of the eq.[19]can be explicitly calculated: 



K, 



/o'-"-(§)' 
n - 1 



A^ In 2 



(22) 



As 1 - n < 0, in grazing incidence the (2//1)'"" term can be 
neglected with respect to f^^"- By isolating the frequency /o and 
using the eq. [20] to derive H{A), we obtain after some algebra, 
for the scattering term of the HEW, 



H{A) = 2 



\6n^K„ 



(«- l)ln2 



sin 6*1 



(23) 



This equation states that: 



1. The H{A) function for a power-law PSD has a power-law 
dependence on the photon energy E oc IjA, i.e., H{E) oc . 
The power-law index y is related to the PSD power-law index 
n through the simple equation: 



r 



n- 1 



(24) 



As 1 < n < 3, y is positive, i.e. H is an increasing function 
of the photon energy. For a fixed value of K„, the HEW di- 
verges quickly for n ^ \ but very slowly for « ^ 3: a PSD 
power-law index close to 2-3 would hence be preferable in 
order to reduce the degradation of focusing performances for 
increasing energies. 

H{A) depends on the sine of the incidence angle at the 
power. In other words, the HEW depends only on the ratio 
sin0i//l: this scaling relation shows that for a given power- 
law PSD (with n < 3) at a given photon wavelength A we can 
reduce the HEW by decreasing the incidence angle. 
H{A) increases with the PSD normalization Kn, as expected: 
the dependence is also a power law with spectral index 



1 



«- 1 



(25) 



As for y(«), the closeness of n to the maximum allowed value 
for fractal surfaces makes less severe the roughness eff'ect on 
imaging degradation. 

The functions yS and y are plotted in fig. [2] For instance, if 
n - 2, y - fi - \, and H(E) increases linearly with both pho- 
ton energy and Kn coefficient. The divergence of indexes p, y for 
« a! 1 makes apparent the importance of obtaining steep PSDs 
in the optical polishing of X-ray mirrors. Finally, it is worth not- 
ing that for n > 3 there is the theoretical possibility of a slight 
decrease of H{E) for increasing energy because y{n) becomes 
negative. 

To clarify the dependence of the HEW on the power-law in- 
dex n and the incidence angle, we depict in fig.[3]and[4]some ex- 
amples of H{E) simulations (single reflection) for some power- 
law PSDs in the photon energy range 0.1-50 keV. The H{E) 
curves were computed using the eq. [23] In fig. [3] the incidence 
angle 9i is fixed at 0.5 deg and the index n is variable; a constant 
n = 1 .8 and a variable 9i is instead assumed in the simulations of 
fig. [4] Note in fig.[3]the slower H{E) increase for larger n and the 
common intersection point, determined by the particular choice 
of the incidence angle and the cr = 4 A value in the window of 
spatial wavelengths [100 H- 0.01 yum]. 
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Fig. 3. li{E) simulations assuming power-law PSDs with con- 
stant cr = 4 A in the spatial wavelengths range [100 H- 0.01//ot], 
but variable power-law index n. The incidence angle is fixed at 
= 0.5 deg. 
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Fig. 4. H(E) simulations assuming a power-law PSD with 
power-law index n = 1.8 and with cr = 4 A in the spatial wave- 
lengths range [100 H- 0.01 //m], but variable incidence angle Q\. 



5. Numerical integration of the PSD 

A power-law PSD is a modelization that can be used for 
optically-polished surfaces. If the polishing process is not op- 
timized or a reflecting layer is grown onto a optically polished 
substrate, several deviations from a power-law trend can be ob- 
served. A typical "bump", for instance, can be present in the PSD 
of multilayer coatings, often in the range of spatial wavelengths 
[10 H- 0.1 yum], as a result of the replication of the substrate to- 
pography and of fluctuations intrinsically related to the random 
deposition process (Spifler et al. 119931 Stearns et al. |1998l l. If 
the PSD deviates significantly from a power-law, the eq.|23]can- 
not be used. However, if the surface PSD has been extensively 
measured over a wide range of spatial frequencies [/m, /m] (wide 
enough to have /m < /()(/!) for all /I), the HEW scattering term 
H{X) can be computed by numerical integration (eqs.[T9landl20b. 
on condition that the following approximation is valid: 



r'p(/)d/^ r 



Pif)df. 



(26) 



The condition above is usually satisfied when /o <K /m i.e. when 
the following inequality holds: 



HiA) <K 



2/t/M 

sin^i 



(27) 



As we are also interested in computing H(A) in hard X-rays 
(small A), there is the possibility that the two integrals in the 
eq.|26] differ by a significant factor. In this case the integral can 
be coiTected by adding the remaining term 



f ' Pif)df 



P(f)df+ P(f)df, (28) 

that can be evaluated, in principle, by measuring the mirror re- 
flectivity within an angular acceptance corresponding to the spa- 
tial frequency /m, and using the Debye-Waller formula to deiive 
cr^; then, the importance of measuring the PSD in a very wide 
frequencies interval becomes apparent. The value of /o depends 
strongly on both incidence angle and photon energy: for soft X- 
rays (< 10 keV) and very small angles (< 0.2 deg) the character- 
istic spatial wavelength I - 1 //o often falls in the millimeter or 
centimeter range. 

It should be noted that, if the detector is small, a fraction of 
the scattered photons can be lost; to account for the finite an- 
gular radius of the detector d (as seen from the optic principal 
plane), one should integrate the PSD over the smaller interval 
[/o, dsin6i/A] to recover the measured H{A) trend. As an alter- 
native method, one can compare the theoretical predictions of 
eg s . [T9] and l20l with the experimental H(A) values, as calculated 
from the Encircled Energy normalized to the photon count fore- 
seen by the Fresnel equations (i.e. with zero roughness), rather 
than to the maximum of the measured Encircled Energy func- 
tion. 



6. Computation of the PSD from the H(A) trend 

If the approach described above can be used to simulate the 
HEW trend from a measured surface PSD, the reverse problem, 
i.e. the derivation of surface PSD from the measured HEW trend 
is also possible. This requires that the figure error contribution 
had been reliably measured, in order to isolate the scattering 
term function H{A) using the eq.|6l 

This problem is interesting for three reasons at least: 

1 . it is a quick, non-destructive surface characterization method 
in terms of its PSD. 

2. The measurement is extended to a large portion of the illu- 
minated optic, hence local surface features are averaged and 
ruled out from the PSD. 

3. For a given HEW(/l) requirement in the telescope sensitiv- 
ity energy band, it allows establishing the maximum allowed 
PSD. 

In order to find an analytical expression for the PSD, we note 
that the spatial frequency /o that scatters at an angular distance 
Hjl from the specular beam is a function only of A, along with 
the eq.|20] Solving for/o, we have 



UA) « H(A)- 



sinft 



2A 



(29) 



We suppose that all scattered photons are collected, so we can 
assume the eq.[T9]as valid. By deriving both sides of eq. [19] with 
respect to A, we have 



d 

dl 



I' 



P(f)df 



In 2 



87r2 sin^ Oi 



A, 



that is. 



In 2 



-|^(l)-^^^^"^=8;r2sin2 0. 



A, 



(30) 



(31) 
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and, using the eq.|29]to compute the derivative of /q: 
-I^(l) + /o^(/o) sm0,dH(A) 



2A dA 



-Pifo) 



In 2 



87r2 sin^ Oi 



(32) 



Now remember that, in grazing incidence, /o <«; 2/i"' by several 
orders of magnitude. Even if P(f) is not a power-law, it is always 
a steeply decreasing function of /. Moreover, it should have over 
[fo, 2A'^] an average PSD index n > 1, for the reasons explained 
in the sect. |4] This means that 



f(/o) 
P(2A-^) 



2A-' \" 2A-' 



fo 



fo 



(33) 



therefore, in practical cases the 2/i 'P(2/l ') term in the eq.[32lis 
negligible with respect to foPifo). Consequently, we can neglect 
the first term of eq. [31] then we have 



d/o 



In 2 



87r2 sin^ Oi 



-A. 



(34) 

]and collecting the constants, we 



Pifo) d 1 


'HiA) 


A dA\ 





0. 



(35) 



Combining this with the eq. |2S 
obtain the final result 

/ 47r2 sin^ 0i 

The eq.^5\enables the computation of the PSD (at the spatial 
frequency given by the eq. l29T l along with the derivative of the 
ratio H(A)I A with respect to A. 

The obtained equation shows that P{f) is inversely propor- 
tional to the derivative of H{A)/A. This result seems strange at 
first glance, because by decreasing H(A) one would obtain a 
larger P(f) (a rougher surface). One should remember, indeed, 
that by reducing H{A) we increase Pifo), but /o is shifted to- 
wards the low frequencies domain, where Pifo) is expected to 
be higher. In fact, the "rough" or "smooth" feature of the surface 
depends on whether fo or Pifo) varies more rapidly, i.e. on the 
overall HiA) trend. 

We can also check the correctness of the eq.l35]bv computing 
the PSD for the particular case of the HiA) derived from the inte- 
gration of a power-law PSD (the eq. |23] derived under the same 
approximation, the eq.[33]l. If the results are correct, the substi- 
tution of the HEW trend of the eq.|23]in the eq.|35]should return 
the original PSD (eq.[2TI). The straightforward, but lengthy cal- 
culation (carried out in appendix |B]) shows that the substitution 
returns 



^(/o)=§ 
Jo 



(36) 



as expected. 

The eq. [35] should be approximately valid also for graded 
multilayers with a slowly-decreasing reflectivity (see sect. 13.2b . 
however, due to the approximations needed to extend the eq. [19] 
to the multilayers, the resulting PSD should be considered a 
"first guess" in this case. Then, the matching of the PSD to the 
required HEW trend should be checked by means of a detailed 
computation of flie XRS PSF(A). 

7. Extension to X-ray mirrors with multiple 
reflections 

The formalism exposed in the previous sections can be extended 
to a double-reflection optic (like a Wolter-I one). In this op- 
tical configuration, photons are firstly reflected by a parabolic 



surface and subsequently by a hyperbolic one. If the smooth- 
surface condition is satisfied, multiple scattering is often negli- 
gible (Willingale 1988) and the scattering diagrams of the two 
reflecting surfaces can be simply summed (De Korte et al. ll98U 
Stearns et al. ,1998 ). The source is assumed to be at infinite dis- 
tance, then X-rays impinge on the two surfaces at the same an- 
gle 6i. If the surface PSDs are the same for both reflections, the 
scattering diagram will be simply doubled. Thus, the integrated 
scattered intensity is also doubled: 



Is = 2IoR] 



1 - exp - 



\6n cr^ sin" 6[ 
^2 



(37) 



The /?F factor is squared in the eq. [37] because each ray is re- 
flected twice: in absence of scattering the reflected power would 
be /o^p, so the half-power scattering angle condition reads 



\6,-6,\ > 



HiA) 



1 2 



and, combining the eqs.l37landl38] we obtain 
167rV2sin2 0i\ 3 



^'^Pl- A^ I 4- 
Solving for cr^, and using the eq.[T6l 

J fa 167r2sin^i 



(38) 



(39) 



(40) 



that differs from the eq. [T9]only in the factor ln(4/3) instead of 
In 2 on right-hand side. Consequently, the corresponding differ- 
ential equation is 



Pifo) d (HiA)\ ^ ln(4/3) 
A dA\ A / 47r2 sin^ 0i 



(41) 



Similar equations can be derived for an optical system with 
an arbitrary number of reflections A^: to compute the HiA) from 
the PSD, 



P(f)df 



2N 



167^2 sin^ ft \2A^-1 



■In 



(42) 



If the PSD is a power-law P(/) - K„/f" we can generalize the 
eq.[23] 



HiA) = 2 



In 



2N 
2N-1 



in- I) 



sin 6\ 

A 



(43) 



note the divergence of the logaritmic factor for increasing A^, 
due to the negative exponent 1/(1 - n). This indicates that HiA) 
increases rapidly with the number of reflections, as expected. 

Finally, we can also generalize the differential eq. [35] to an 
arbitrary number of reflections. 



Pifo) dlHiA)\ ln(#T) 



d lHiA)\ 
dA\ A j 



0. 



A dA\ A } 47r2sin^fti 
In the eqs.l42]andl44l fo is always related to HiX) by the eq 



(44) 
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Fig. 5. an hypothetical PSD with reasonable values and a PSD 
break around a 100 jum spatial wavelength (dashed line). This 
PSD is adopted to compute the corresponding HEW trends for 
1,2,3 reflections at a 0.3 deg grazing incidence angle (fig.©. The 
achieved HEW trends were used to re-calculate the respective 
PSDs (solid line). For clarity, we do not plot the single PSDs, 
but just their overlap. 



8. An example 

As an application of the equations reported above, we shall 
make use of a simulated surface PSD with reasonable values, 
that is not a power-law. The PSD (see fig. |5] dashed line) is ex- 
tended from 10^ fim down to a 0.01 ;um spatial wavelength, with 
a break around lOOyum: at the lowest frequencies the PSD is 
steep (n ^ 2.3), whereas at the largest frequencies it is smoother 
(n =s 1.3). From the discussion in sect. |4]concerning the relation 
between the exponents of the PSD and the HEW (eq. l24l l. we 
should expect that the PSD break causes a slope change in the 
function //(/I): however, as the actual PSD is not a power-law, 
the H(A) function should be computed by means of the eqs. |20] 
and|42] Before carrying out the integration, we can remark qual- 
itatively that, as we increase the photon energy, the highest fre- 
quencies in the PSD (where the PSD index becomes smaller) 
become important; hence, we can expect a steeper increase of 
the HEW at the highest energies. 

The analysis is made quantitative in fig. |6l where we show 
the calculated HEW trends from the PSD in fig. |5](the dashed 
line) by means of the eqs.|20]and|42] assuming 1,2,3 reflections 
at the same grazing incidence angle (0.3 deg). The approxima- 
tion of eq.|26]was adopted. In addition to the scattering term, 15 
arcsec of HEW due to figure errors were added in quadrature. 
The HEW increases slowly (concave downwards) at low ener- 
gies, corresponding to a frequency /o in the steeper part of the 
PSD. Then it increases more steeply (concave upwards) when 
the energy becomes large enough to set /o in the portion of the 
spectrum with n ^ 1.3. By increasing the number of reflections, 
the HEW values also increase, and the "turning point" where the 
HEW starts to diverge (arrows in fig. ^ shifts at lower X-ray 
energies. All the calculation is based on the assumption that the 
contribution of the PSD over the maximum measured frequency 
/m - 0.01 jum is negligible. Otherwise, the computed HEW val- 
ues will be underestimated (see sect.|5l). 

In addition to the general trend of the HEW, there are oscil- 
lations due to small irregularities in the adopted PSD: the cal- 
culation is, in fact, very sensitive to small variations of the PSD 
values. Notice that for a definite energy all the frequencies larger 
than /o contribute to the HEW value, even if the largest contri- 




5 10 15 20 25 30 35 40 45 50 
Photon Energy (keV) 

Fig. 6. the HEW trend computed from the PSD for 1,2,3 reflec- 
tions, plus 15 arcsec of HEW due to figure errors. The HEW 
trends were used to compute back the PSD (the solid line in 
fig E) to verify the reversibility of the calculation. The energy 
at which the concavity change takes place is also indicated (ar- 
rows). 

bution comes from frequencies near /q: this is a consequence of 
the steeply decreasing trend of the PSD. 

We checked the reversibility of the result by computing the 
PSD from the HEW trends (after subtracting in quadrature 15 
arcsec figure error) by means of the eq. |44] with the respec- 
tive value of A^. The resulting PSDs (the solid line in fig. |5]) 
were overplotted to the initial PSD, with a perfect superposition. 
Each obtained PSD has, indeed, an extent of spatial frequencies 
smaller than the initial one: the overall PSD ranges from 10^ 
to 11 yum (vs. the initial 10^ 0.01 yum), and the smaller wave- 
lengths could be computed from the HEW trend with N = 3. The 
limitation in spatial frequency ranges occurs for two reasons: 

1. small f - large i. all the power scattered by the lowest fre- 
quencies is found at angles less than 1/2 HEW even for the 
lowest energies being considered: therefore, that part of the 
spectrum is not necessary to compute the HEW in the energy 
range of interest; 

2. large f - small I: the PSD is computed from a derivative, 
therefore the information concerning the absolute magnitude 
of the HEW is substantially lost. This information was in- 
cluded in the integral of the PSD (eq. l42l l for the maximum 
considered energy. 

Therefore, from the integral in the eq. |42] we cannot re- 
cover the PSD over the minimum computed spatial wavelength 
(llyLfm, using the HEW trend with A^ = 3), but we can at 
least calculate the value of cr at spatial wavelengths smaller than 
1 1 /urn. Substituting the incidence angle and the minimum pho- 
ton wavelength being considered (A - 0.24 A) in the eq. |42] 
with A = 3 and with the approximation of the eq. |26] we ob- 
tain cr = 1 .6 A, in perfect agreement with the value computed 
from the original PSD. 

Summing up, for a given incidence angle the H(A) function 
in a definite photon energy range is equivalent to the PSD in 
a coiTesponding range of spatial frequencies / (or equivalently, 
spatial wavelengths I), plus the integral of the PSD beyond the 
maximum frequency being computed. Therefore, requirements 
of a definite HEW(/l) function in designing an X-ray optical sys- 
tem can be translated in terms of PSD in a frequencies range 
[/min, /max] plus the surfacc rms at frequencies beyond /max- The 
usefulness of such a relationship is apparent. 
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9. Conclusions 

In the previous pages we have developed useful equations to 
compute the contribution of the X-ray scattering to the HEW of 
a grazing incidence X-ray optic, by means of a simple integra- 
tion. The formalism has been inverted in order to derive the PSD 
of the surface from the function H(A), and it can be extended 
to an arbitrary number of reflections at the same incidence an- 
gle. The equations are valid for a single-layer coating miiTor, but 
they can be approximately applied to a multilayer-coated mir- 
ror. This approach is particularly useful in order to establish the 
surface finishing level needed to keep the X-ray scattering HEW 
of X-ray optics within the limits fixed by the X-ray telescope 
requirements. 

It should be remarked that the reasoning was developed for 
the Half-Energy Width, but it can be extended to any angular 
diameter including a fraction rj of the energy spread around the 
focal point. To do this, it is sufficient to substitute the logarithmic 
factors in equations l42l |43l l44l 



In 



2N 
\2N- 1 



In 



A^- 1 +77 



(45) 



and for instance, to compute the 90%-energy diameter for a dou- 
ble reflection mirror, simply substitute rj = 0.9 and N - 2. The 
proof is straightforward: however, one should always keep in 
mind that the energy diameters computed with this method can 
be considered valid only if they are much smaller than the inci- 
dence angle 0i. 

Notice that in the development of the exposed formalism we 
have supposed, in addition to the smooth-surface condition, two 
additional hypotheses: 

1 . the source is at infinite distance from the mirror 

2. the X-ray detector is large enough to collect all the scattered 
photons. 

In order to apply the mentioned equations to experimental 
calibrations of X-ray optics at existing facilities (like MPE- 
PANTER), where the source is at a finite distance and the de- 
tector has a finite size, some corrections should be taken into 
account. We will deal with their quantification in a subsequent 
paper. 



Appendix A: Extension to multilayer coatings 

Here we provide with a plausibility argument to extend the for- 
malism of sect. l3.1l to mirror shells with multilayer coatings (see 
sect. 13.21 ). The intensity of a scattered wave at each interface is 
proportional to its PSD as per the eq.[TT] and the overall scatter- 
ing diagram will be their coherent interference. To simplify the 
notation, we neglect the X-ray refraction and we suppose that 
the incidence angle is beyond the critical angles of the multilayer 
components. The electric field scattered by the A:* interface can 
be written as 



Ek = EQTkri,Xk(f)exp(-i(pk), 



(A.l) 



where is the single-boundary amplitude reflectivity, Eq the 
incident electric field amplitude, the weights Tk are the rel- 
ative amplitudes of the electric field in the stack (in scalar, 
single-scattering approximation), and account for the extinction 
of the incident X-rays due to gradual reflection and absorption. 
Xk(f) is the single-boundary scattering power (proportional to 



the PS D(f) amplitude), and (p^ is the phase of the scattered wave 
at 0s by the A;* interface 



sin 6i + sin 6s 
h - 2n ■ Zk, 

A 



(A.2) 



where Zk is the depth of the A:* interface with respect to the outer 
surface of the multilayer. Now, the measured intensity is 



\Esca 



I^Ek 

k=Q 



m^iXkiDi" 



k=o 



(A.3) 



Now, \E()\^ - /(), the incident X-ray flux intensity, and |Xj;(/)p is 
proportional to the interfacial PSD P{f), which is independent 
of k by hypothesis. Assuming the proportionality factor of eq.fTTI 
for |Xt(/)p, we obtain for the scattering diagram 



1 dh 



and if we set 



sin^ OiPif) 



k=0 



(A.4) 



k=() 



(A.5) 



the eq. IA.4I becomes analogous to the eq. (TT] with K^{Oi,0s) 
playing the role of Rp. Note that K,i(6i,6i) - RA(6i), the mul- 
tilayer reflectivity in single reflection approximation. As before, 
we write the scattering diagram for a mirror with axial symmetry 
as a function of the angular distance from the focus a = \0i - 0^1 
averaging the contributions of negative and positive frequencies 

T^^%r sin' 0iPif)[KA(0i, Oi-a) + KM, 0i + a)] . (A.6) 
i() at/s 

For a single reflection optic, we can calculate the scattered power 
over H/2, where H is the scattering term of optic Half-Energy 
Width: 



h[a > H/2] = ^loRAiOi). 
Now, the steps [TS] and [T9l can be repeated: 
' KM, Oi + a) + KMi, Oi-a) _ 



RMi) 



-Pif)df- 



In 2 

8;r2 sin^ 0; 



(A.7) 



(A.8) 



where /o is still defined by the eq.|20] For small scattering angles 
(a <sc Oi), since we assumed a slow variation of R i over angular 
scales of H/2 (and the same occurs for K^), we can approximate 



KA(6i,6i±a) ^ RA{0i)±a 



dKAi0i,0s) 



d0s 

Substituting in the eq. IA.8I we obtain 
A^\n2 



r 



P(f)df^ 



16n^ sin^ ft 



(A.9) 



(A. 10) 



because the two derivatives have opposite sign and cancel out. 
This is the same equation found for the case of a single-layer 
coating (eq.[T9]). 
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Appendix B: Derivation of the PSD from the IHEW 
for a fractal surface (single reflection) 

We recall here the H{A) trend for a power-law PSD (eg. [23]): 



HiA) = 2 



167T^K„ 



(n- l)ln2 









fsinflij 



(B.l) 



we verify that it returns a power-law PSD if substituted in the 
differential eq. [35l 



P(fo) d lHiA)\ ^ In 2 ^ ^ 
A dA\ A I 47r2sin3 6li 



(B.2) 



To simplify the notation, we write simply H instead of H(A): by 
carrying out the derivation, 



1 d ///\ _ 4 
A dAXj) ~ ~ n-\ 



in - l)ln2 
Using again the eq. lB.ll 

1 d lH\ 2H , 



(sinfl,)^^^ 



--{-) 

A dA\Al 



n-l 

hence, the related PSD is 



P(fo) = - 



In 2 



47r2 sin^ 0i 



A dA\Al 



A^\n2 n-l 



Now, we can derive (n - l)/2 from the eq. lB.ll 
n - 1 AyrHKn /Hy" Isin 6/ ^ 



2 ln2 \2/ \ ^ 

and combining the eqs. lB.5lfB.6l one obtains 



P(fo) = K„ 



2A 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



(B.7) 



(B.8) 



that is, by recalling the eq. |29] 

P(fo) = § 
Jo 

i.e., the expected power-law PSD. 
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